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Abstract. The statistics of nested spiral self— avoiding loops, which is closely 
related to the partition of integers into decreasing parts, has been studied on the 
square and triangular lattices. The number of configurations with N steps is cjv — 
(y2/24)Af-3/2 exp(7ry^A^i/2) ^^^^^ ^j^gj^. average size Xn ^ {l/2n)-y/^N^/^ IniV to 
leading order on the square lattice while the corresponding values for the triangular 
lattice are cn (33/4/l6)Af-5/4 exp((7r/\/3)Afi/2) and Xn ~ l/(7rVS)N^/'^ \nN. 



Some years ago, the number of A^-step spiral self-avoiding loops have been 
calculated for the square and triangular lattices (Manna 1985, Lin et al 1986). These 
works followed the introduction of the spiral self-avoiding walk (Privman 1983) for 
which a lot of exact results were obtained by a succession of authors (Blote and Hilhorst 
1984, Whittington 1984, Gutmann and Wormald 1984, Joyce 1984, Guttmann and 
Hirschhorn 1984, Lin 1985, Joyce and Brak 1985, Lin and Liu 1986). While the 
number of spiral self-avoiding loops grows with like a nonuniversal, i.e. lattice- 
dependent power, the number of spiral self-avoiding walks also behaves in an unusual 
way 

Cj^ ~ AN-^ exp (AiV^/2) (1) 

where both 9 and A are lattice-dependent quantities. It follows that the asymptotic 
entropy per step decays as N~^^'^ instead of giving a constant like in the ordinary 
or directed self-avoiding walks. Following the work of Privman, a close connection 
between this problem and the theory of partitions of integers (Andrews 1976) was 
noticed (Derrida and Nadal 1984, Redner and de Arcangelis 1984, Klein et al 1984). 

In this letter we present some exact results concerning the statistics of nested 
spiral loops which are self- and mutually avoiding and piled up around a site chosen 
as origin. We study such spirals on the square lattice where on a loop only 90° turns 
in the same direction are allowed so that the loops are rectangular-shaped (figure 1) 
and on the triangular lattice where the restriction to 120° turns leads to triangles 
among which one only keeps those pointing up (figure 2). 
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Figure 1. Nested spiral self-avoiding loops on the square lattice: with 90° turns in 
the same direction, rectangular-shaped loops are obtained. Each step is assigned a 
weight z and the size is ^ = ^2k—i nested-loop configuration corresponds 

to four independent partitions of integers into decreasing parts numbered 1 to 4 and 
each partition is duplicated (heavy lines). 




Figure 2. Nested spiral self-avoiding loops on the triangular lattice: with 120° turns 
in the same direction, triangular-shaped loops are obtained. Each step is assigned a 
weight z and the size is ^ = "^^—1 ^k- The nested-loop configuration corresponds 
to three independent partitions of integers into decreasing parts numbered 1 to 3 and 
each partition is triplicated (heavy lines). 



Let US introduce the generating function 

G^(z,a;) = 5^c^(L,X)A-^ (2) 

N,X 

for the number of configurations with N steps, L loops and size X — '^^^-^Xj^ where 
X is the distance from the origin to the Lth loop. On the square lattice, with the 
notations of figure 1, one may write 

00 00 xi—l 00 00 X2—X1 — I 

xi=2yi=2 Xi = l X2=xi+2 y2=yi+'2 X2=l 

00 00 xl—^l-i — I 

■■■ E E (vL-yL-i-i) E ^'^-^^^'=+^'=^6'^^.^^'''= (3) 

xl=xl-i+2 yL=yL-i+2 Xl=1 
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Introducing 



^k = yk- Vk-i 2/0 = 



one obtains 



L 



(4a) 
(4b) 



(5) 



k=l k=l 

and the generating function may be factorized as 

L 



k=l 



Z (^4^) EK-i) 



,2(L-fc+l)nfe 



(6) 



or 



G1 



)3(1 -^^fcgc) 



(7) 



In the same way, on the triangular lattice with the notations of figure 2, one has 



CO x\—2 00 X2—X1—2 

Gi^(^,a;)= J] J](xi-Xi-1) Yl E (^2-^i-^2-l)--- 

a;i=3Xi = l X2=xi+3 X2 = l 

CXD xl—xl-1—2 

... J2 E (^L-^L-i-^L-l)^'^-^"'=e-S.^i^M8) 



SO that using (4a) 



k=l 



.mfe=3 



^ _ 2)e- - (m, - l)e^- + e--'»'= ..3(L-fe+i)m.' 

(1-6^^)2 



(9) 



and finally 



(10) 



These results may be obtained more directly by noticing the connection with the 
number of partitions of integers into L decreasing parts which is illustrated on figure 1 
for the square lattice and figure 2 for the triangular lattice. Nested loops configurations 
are in one-to-one correspondence with four independent partitions in the first case 
and three independent partitions in the second case. The generating functions for the 
ncstcd-loop problem are directly obtained as powers of the generating function for the 
partition of integers into L decreasing parts (Andrews 1976) 



g^it) ={t + i' + + ...){£' + 1"" + t'' + ...).. .{t^^ + 1''^ +1^"^ + .) 



^ 4.k 

fe=i 



(11) 



t^ 
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multiplied by a modified generating function for which a partition of size X is, as in 
equations (3) and (8), weighted by e'^-^ so that 

g^it, ij) ^{te^ + t^e^^ + ■ ■ ■){t^e^ + t^e""^ + ■■■)■■■ {t^e^^ + t^^e^'" + ■ ■ ■) 

= rr ^^^^ 

-I- -I- 1 - e'^t^ 

k=l 

One may verify that the mapping requires t = for loops on the square lattice 
(figure 1) and t = z^ on the triangular lattice (figure 2); it follows that 

G^{t,uj) = [gj^{t)r^ g^{t,u)\,^^, (13) 

with p = 4, g = 2, on the square lattice and p = 3, g = 3, on the triangular lattice in 
agreement with equations (7) and (10). 

Let us first consider the generating function G{t) for the number of configurations 
CjY = 5^^ ^ Cj^f{L,X) with N steps and any number of loops which is given by 

oo oo 

Git) = Gl(^,^)l=o = E ^9Lit)r it=.. (14) 

L=l L=l 

The bchaviolu' of for large N values is governed by the behaviour of ^^^(t) in the 
vicinity of its singularity at t = l. With t = l — r] and r]^0^, the main contribution 
to the sum in equation (14) comes from values of L near LQ = ln2/?7 for which gi^if]) 
is maximum in L (Blote and Hilhorst 1984, des Cloizeaux and Jannink 1987) and an 
expansion of \ngj^{rj) near Lq (see the appendix) leads to 

ln<,.W-^-iln(^^)-i(L,-ln2f (18) 

SO that the sum over L in G{t) may be transformed into a Gaussian integral and one 
obtains 

oo p-1 / 2 \ 

Git) = E expbln^^(t)] ^ (2p)-^/2 (^) ' exp (^^J (16) 

for the leading contribution. With n = N/q, the generating function may be written 
as 

oo 

G{t) = E ^^=.n^" (17) 
n=0 

SO that the number of configurations is given by the Cauchy formula 

1 / .G{t) 



where (C) is a circle of radius r < 1 centered at the origin. The integral may 
be evaluated using the saddle-point method by deforming the contour. Through a 
Gaussian integration near the saddle-point at t = l— 7r(p/12n)^/^ one obtains 



c 



N 



2P+I3P/4 



Ar-4^exp(7rjf iV^/M (19) 
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to leading order. With the appropriate p and q values one obtains 



24 \ V 3 



^^-^Ar-3/2exp(7r^/jArV2| (20a) 



33/4 



16^- ^"V^ 



Ar-5/4exp(^iVV2) (20b) 



For both lattices, the exponential term is the same as for the outward spiral self- 
avoiding walk (Blote and Hilhorst 1984, Joyce 1984, Guttmann and Wormald 1984, 
Joyce and Brak 1985, Lin and Liu 1986) whereas the power of N in the prefactor is 

different. 

According to equation (2), the mean size for L-loop configurations is given by 



and may be rewritten as 



ui=0 



fc=i 



k=Q ^ ^ 

where the first term in the sum should be understood as the limit of the ratio when 
A;— >0. When L is unrestricted, the mean size becomes 

Changing the sum over L into an integral, the value LQ = ln2/?7 corresponding to the 
maximum in L of Gj;^{t) is selected and 

°° j-fc/1 — k\ 
A;=0 

Putting apart the first term and rearranging the sum, one obtains 



(25) 



^, , ln2 ^ / 1 1 \ 

The sum may be evaluated using the Euler-Maclaurin formula and 

X{rj)c^-H-) (26) 
V V 



to leading order. On the other hand the number of steps reads 



dlnG{z'^) _ pgTT^ 



dliiz 



zi=l-r) 



(27) 
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so that 



X ~ 



1 /T2 
2tt V M 



(28) 



and 



1 . 3„i/ 



7rv3 



(29a) 
(29b) 



One recovers the characteristic behaviour of the spiral self-avoiding walk (Blote and 
Hilhorst 1984, Liu and Lin 1985) with an exponent v = 1/2 and a logarithmic 
correction. 



Appendix 

Using equation (11) with t = l—r] in the limit 77^0+ one has 



In^^(t) = - In it-'' - 1) ~ - ^ [In {e''^ - 1) - In/cr/] -J^^nkri 



k=l k=l 

The first sum may be replaced by an integral using 



k=l 



V J n. 
' 2 



24 



1, 



SO that 



1 r(L+^)v 

lngM = / [ln(e" - 1) -lnw]dw-Llnr/-ln(L!) + 0(r/) 

V Jo 



(Al) 



(A2) 



(A3) 



Integrating the second term and using Stirling formula, one gets to leading order 



1 1 /2tt\ 

IngM ^ — / In (e" - l)du - In - 

V Jo ^ V ^ / 



(A4) 



Expanding the remaining integral, considered as a function of its upper limit, near 
Uo = -^o^~^n2 with 



/■In 2 

/ In (e" - 
Jo 



l)du = - 



12 



(A5) 



one finally obtains 



ln^^(r7)^^-^ln(^)-^(L,y-ln2)2 



(A6) 
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